PROCEEDINGS 

OF SCIENCE 



A Ghost Story: Ghosts and Gluons in the IR regime 
of QCD 



o 
o 

(N 
> 

o 



Or 
43 



> 
00 

O 



On 
O 



X 



Olivier Pene* 



E-mail: olivier .pene@th.u-psud. fr 



Ph. Boucaud 

LfT Or say (CNRS) 

E-mail: |philippe . boucaud@th . u-psud . f r 



J. P. Leroy 

LPT Orsay (CNRS) 



E-mail: jean-pierre . leroy@th . u-psud. f r 



A. Le Yaouanc 

LPT Orsay (CNRS) 



E-mail: alain . le-yaouanc@th . u-psud . f r 



J. Micheli 

LPT Orsay (CNRS) 



E-mail: jacques . micheli@th . u-psud . f r 



J. Rodriguez-Quintero 

Dpto. Fisica Aplicada, Huelva $ 
E-mail: 



jose . rodriguez@df aie . uhu . es 



We discuss the different methods to obtain reliable informations about the deep infra-red be- 
haviour of the gluon and ghost Green functions in QCD. We argue that a clever combination of 
analytical inputs and numerical ones is necessary. We illustrate this statement about the distinction 
between two classes of solutions of the ghost propagator Dyson-Schwinger equation (GPDSE). 
We conclude that the solution II ("decoupling") with a finite renormalised ghost dressing function 
at zero momentum is strongly favored by lattice QCD, We derive a method to solve numerically 
the GPDSE using lattice inputs concerning the gluon propagator. We derive an analytical small 
momentum expansion of the Ghost dressing function. We prove from the large cut-off behaviour 
of the ghost propagator renormalisation constant, Z3, that the bare ghost dressing function is infi- 
nite at the infinite cut-off limit. 

International Workshop on QCD Green 's Functions, Confinement and Phenomenology 
September 7-11 , 2009 
ECT Trento, Italy 



* Speaker. 

' Laboratoire de Physique Theorique, Unite Mixte de Recherche 8627 du Centre National de la Recherche Scien- 
tifique Universite de Paris XI, Batiment 210, 91405 Orsay Cedex, France 

''Fac. Ciencias Experimentales, Universidad de Huelva, 21071 Huelva, Spain 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribntion-NonCommercial-ShareAlike Licence. 



http: //pos . sissa . it/ 



A Ghost Story: Ghosts and Gluons in the IR regime of QCD 



Olivier Pene 



UHU-FP/08-010 
LPT-Orsay/09-88 

1. Introduction 

During the ten years ranging from 1965 to 1974 was invented the standard model of particle 
physics. This was a major and often overlooked scientific event. During the early sixties, the four 
fundamental interactions where known, but the weak interaction was only described by Fermi's ef- 
fective theory and the strong interaction seemed to be even further away from any sound theoretical 
description, precisely because, being strong, it seemed impossible to control by expanding around 
a small parameter, as had been the case for quantum electrodynamics. 

Then came the miracle. A quantum field theory containing quarks and gauge particles named 
gluons was proposed and its major property was isolated: asymptotic freedom. It was a miracle 
because its formulation is extremely compact, with only rif + 1 free parameters 1 naming n f the 
number of quark flavors, i.e. tif = 6, the beautiful constraints of gauge symmetry, while the field 
of its applications is huge. Up to now no strong argument has been presented which could allow to 
deny QCD to be the theory of strong interactions. Of course there are drawbacks. The first is that 
we are not able to extract very accurate predictions from QCD's premises. 

But the most frustrating unsolved problem is the inexistence of a real proof of the confinement 
property, i.e. of the observation that only hadrons are observed in nature and never isolated quarks 
or gluons. We are all convinced that confinement is a property of QCD. Confinement is an experi- 
mental fact. Furthermore lattice-QCD (LQCD) calculations, which are based on QCD's principles, 
provide results in full agreement with confinement. But this is not a proof. 

Confinement is the major issue of this meeting and we all believe that it has to be looked 
for in the infrared behaviour of QCD. We will hear in this conference discussions around several 
approaches to confinement. Our approach will not be to follow or criticise some confinement sce- 
narios, but rather to try to provide reliable answers to the question: How do Green functions behave 
in the deep infrared. In this talk we will, for the sake of simplicity, restrict ourselves to the quark- 
less pure Yang-Mills theory. We assume that the main features of QCD's infrared properties are 
present in Yang-Mills, at odds with Gribov's hypothesis that the light quark supercritical binding 
was the origin of confinement [fflj. 

1.1 Tools to handle the Green functions in the deep infrared 

There exist analytical tools which are mainly Ward-Slavnov-Taylor identities (WSTI) and 
Dyson-Schwinger equations (DSE). There exists a numerical tool which is LQCD. 

WSTI and DSE are exact. They can be derived rigorously from the path integral formulation 
of QCD. However, WSTI's are a necessary a posteriori check but do not constrain so much while 
DSE's are a very large set of coupled non linear integral equations. Trying to solve the latter is a 
formidable task and it is not clear how many solutions exist. There can even be an infinity of them. 

LQCD is exact, it is really an approximation of QCD, however it is only numerical, leading to 
an intrinsic uncertainty, and, as we already mentioned, the accuracy is poor. 

1 rif + 2 if we count the strong CP term. 
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We believe that it is extremely fruitful to combine both the analytical and numerical ap- 
proaches. Indeed we will see an example in which LQCD allows to decide between two very 
different classes of solutions of DSE's. And next we will see an example in which analytic meth- 
ods provide the functional form of the ghost propagator in the deep infrared, thus allowing to 
extrapolate to zero momentum, where no direct LQCD calculation is possible. 

1.2 Notations and definitions 

In latin languages the translation of "ghost" starts with an "F", while "gluons" starts with a "G" 
in all languages we know. Therefore we use the following notations: the bare gluon propagator is 
written 



n ab / 2 A 2x _ G(p 2 ,A 2 ) 



p 



'ab 



x p v- Pv 

P 



(1.1) 



where G(p 2 ,A 2 ) is the bare gluon dressing function, A is the ultraviolet cut-off, inverse lattice 
spacing a~ l in the lattice case. The bare ghost propagator is written 

P 

where F(p 2 ,A 2 ) is the bare ghost dressing function. 

The corresponding renormalized quantities are labelled by the R subscript: 

G R (p 2 ,li 2 ) = lim Z 3 -V,A) G(p 2 ,A 2 ) F R (p 2 ,^ 2 ) = Urn Z^(^i 2 ,A) F(p 2 ,A 2 ) (1.3) 

where /J. is the renormalisation scale. 

An important remark for the following is that [||, |3|] 



F(0,A) = Z 3 (M 2 ,A) (F R (0,ii 2 ) + 
G(0,A) = Z 3 (At 2 ,A) (^^(O^ 2 ) + 



1 

A 2 " 
1 

A 2 " 



(1.4) 



In the MOM renormalisation scheme, the renormalised quantities are set equal to their tree 
value when the momentum is equal to the renormalisation scale: 

G R ^ 2 ,n 2 )=F R (n 2 ,^ 2 ) = l (1.5) 



whence, using eq. (1.3) 



Z 3 Gu 2 ,A 2 )=G(At 2 ,A 2 ), Z 3 (At 2 ,A 2 )=F( i u 2 ,A 2 ). (1.6) 
The bare ghost-ghost-gluon vertex is parametrised by 

ff c (-q,k;q-k) = ig f abc ( q v Hi(q,k) + (q - k) v H 2 {q,k) ) , (1.7) 
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where k (q) is the incoming (outgoing) ghost momentum. Taylor's theorem [Q] implies that the 
ghost-ghost-gluon vertex becomes trivial when the incoming momentum vanishes 

H l (q,0)+H 2 (q,0) = \. (1.8) 

This implies that if we take this kinematics to renormalise the ghost-ghost-gluon vertex, the vertex 
renormalisation constant is z\ = 1. 

Exploiting this property we can define very simply "Taylor's coupling constant" by ^ 

a T (p 2 ,A 2 ) = F 2 (p 2 ,A 2 )G(p 2 ,A 2 ). (1.9) 

where gl(A 2 ) is the bare coupling constant. Notice that, while gQ,F and G depend logarithmically 
on the cut-off A, a F only depends on it via inverse powers 0(1/ A 2 ). 

Finally we assume, as everybody does, some simple power law in the deep infrared: 

G(p 2 ,A 2 ) ~SP 2 ) aG *V,A 2 ) ~> 2 r a T (p 2 ,A 2 ) ~Ap 2 ) 2aF+aG (1-10) 



2. Two classes of solutions to the ghost propagator Dyson-Scwinger equation 

Let us consider the ghost propagator Dyson-Schwinger equation (GPDSE). It was claimed by 
many authors trying to solve the DSE's that a general conclusion was that 2a F + ccc = or, in other 
words that a F (p 2 ) — > ct > when p 2 — » 0. On the other hand many indications from lattice QCD 
show a strong vanishing of a F (p 2 ), see a recent result at very small momenta in fig. 5 of [gj. 

Looking into details of the GPDSE we found that there were indeed two classes of solutions [{7|] 

• Solution I: 2a F + a G = 0, a T (p 2 ) r^ ct > and a F < 0, F(p 2 ,A 2 ) oo 

• Solution II: a F = 0, F(p 2 ,A 2 ) ct > and, using the lattice evidence that a G > 0, 2a F + 

a G > 0, a T {p 2 ) — > 

This is valid at fixed cut-off A. Similar conclusions hold for the renormalised quantities. Solution 
I is often called the "scaling solution" while solution II is called for some reason the "decoupling 
solution". 

2.1 Schematic proof of the existence of the two solutions 

The GPDSE writes in our notations as 



1 



F{k 2 ,A) 



d A q F(q 2 ,A)G({q-k,A) 2 



(2n) 



q 2 {q-kf 



k 2 



tfi(?,*,A) , (2.1) 



Using Taylor's theorem eq. (IO), completed with indications from perturbative QCD, we assume 



that in eq. (2.1) H\ is regular: never vanishing nor infinite, and not too far from 1. In practice we 



will take it as a constant close to 1. This hypothesis is rather usual. The r.h.s of eq. (2. 1 ) is divergent 
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at fixed A since the integrand GF decreases at large q 2 as a 35 / 44 which is not enough to make the 
integral convergent. Therefore we prefer to regularize it by using a subtracted GPDSE: 



1+glNc 



F(k 2 ) F{k' 2 ) 




(2.2) 



where we have assumed Hi to be a constant and omitted to write A. The r.h.s. now is convergent. 
Let us assume that we rescale all momenta by a common factor A — > 0, we count the power be- 



haviour of the l.h.s and the r.h.s. From eq. ( 1.10 ) the l.h.s behaves like (X 2 )~ aF and the r.h.s as 
(X 2 ) aa+aF . Matching both sides leads to solution I: 2ap + OCg = 0. However, there is a loophole 
in this argument when \/F{p 2 ) — > l/F(0) > 0, i.e. when ctp = since then the l.h.s. vanishes to 
leading order. To get a relation we need to go to the subleading behaviour of 1 /F(p ). We are then 
in the case of solution II: ctp = and no constraint on 2ap + 0Cg- This is the proof. More details 
can be found in 

2.2 Numerical resolution of the GPDSE 

In order to understand better the relationship between these two classes of solutions we have 
performed a numerical solution of the GPDSE Since we consider only one DSE, we need 
additional inputs. Our inputs are: 

• The gluon propagator is taken from lattice QCD. It is extrapolated to the large momenta 
using perturbative QCD formulae and to zero momentum assuming a finite, non zero limit, 
as strongly indicated by lattice QCD. 

• The ghost-ghost gluon vertex is taken to be constant as justified above from Taylor's theorem. 

• The coupling constant multiplies the vertex function which we assume to be constant. This 
product, a rescaled coupling constant, is taken as a free parameter. 

We then fit this parameter to recover a solution in agreement with the ghost propagator computed 
by lattice QCD. This exercize can be performed, mutatis mutandis, with bare quantities or renor- 
malized ones. In the latter case we define the rescaled coupling constant by 

g 2 = N c g\z\H\R (2.3) 

where z{ is the vertex renormalisation constant and our renormalisation scale is chosen to be 1.5 
GeV. 



Our result is that there is one critical value of the rescaled coupling constant g c =33.198 for 
which the renormalised ghost dressing function diverges at zero momentum, solution I ("scaling"), 
while for all smaller g 2 , F(0) is finite, solution II ("decoupling"). Fitting to the values of F(k 2 ) 
from lattice data gives g 2 = 29. The plots are shown in fig. ||. Not surprisingly, the plot ||| shows 



that the product F G, proportional to Taylor's coupling constant, eq. ( |1.9| ), goes to a constant for 



the critical g c and vanishes for any smaller g , fitting lattice data for g = 29. 
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2.5 



dressing function lattice result 
best SD solution 
singular SD solution 




Figure 1: Comparison between the lattice SU(3) data at /3 = 5.8 and with a volume 32 4 for the ghost 
dressing function and our continuum SD prediction renormalised at jj, = 1 .5 GeV for g 2 = 29. (solid line) ; 
the agreement is striking ; also shown is the singular solution which exists only at g 2 = 33.198.... (broken 
line), and which is obviously excluded. 



2.3 Expansion of the ghost propagator at small momentum 

From the GPDSE one can derive the low momentum expansion of the Ghost dressing function 
in the case of solution II [J8|] . If we assume that the gluon propagator goes to a finite constant, which 
implies that ar(k 2 ) °= k 2 , it takes a simple form: 



F R {k 2 ^ 2 ) ~ F s (0 )A t 2 ) (l + ^|^ar(* 2 )log(* 2 ) 
F(k 2 ,A) * F(0,A) (\ + ^a T (k 2 )log(k 2 ] 



(2.4) 



This formula, which can be refined [@], is very useful since it allows an extrapolation of lat- 
tice data down to zero momentum. This is an exemple in which an analytic method supports the 
numerical one. 



3. What do we learn from lattice QCD 

This will be a very brief section as everything has been covered in Teresa Mendes's talk. 
Let us just mention recent publications, which present results obtained with particularly large vol- 
umes and thus small momenta. What follows concerns bare Green functions at some finite cut-off. 
Cucchieri-Mendes have studied the SU (2) case [JlOp: their fig.2 shows a bending of the ghost dress- 
ing function perfectly compatible with solution II ("decoupling"). In [11] they consider the j3 = 
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Figure 2: Comparison between our lattice SU(3) data at j3 = 5.8 and with a volume 32 4 for the product of 
gluon times ghost square dressing functions GR(k)FR(k) 2 , renormalised at jx = 1.5 GeV, and the correspond- 
ing curve for the continuum singular solution (Xq + Iolf = 0, which exists only at g 2 ~ 33.198, obviously 
excluded. Also shown is our continuum regular solution for g 2 = 29 (solid line) for which the agreement is 
striking. 



situation and exhibit bounds on the gluon propagator (their fig. 4). Bogolubsky et al. [g] consider 
SU (3): their fig. 2 shows that the gluon propagator goes to a non zero constant at zero momentum, 
fig. 4 shows also a bending of the ghost dressing function and fig. 5 clearly shows a vanishing of 
a,j at zero momentum. The general conlusion is that the gluon propagator goes to a non zero con- 
stant, the ghost dressing function may go to a finite non zero limit, and Taylor's coupling constant 
clearly vanishes at zero momentum. If the finiteness of the ghost dressing function is today only 
an indication, the vanishing of the coupling constant is compelling, thus contradicting solution I 
("scaling"). Now, since the analytic GPDSE method says that there exists only these two classes 
of solutions, we may conclude that nature has chosen solution II and thus that the ghost dressing 
function must indeed go to a finite non zero constant at finite cut-off. As we see this is an exemple 
in which the LQCD numerical method allows to discriminate between two classes of solutions of 
the GPDSE. 

One remark is in order here. We use the denomination "coupling constant" in a very general 
sense: any well defined quantity which in the ultraviolet is asymptotically equivalent to, say, 0Cj^, 
is eligible for the denomination "coupling constant". (Xt defined in eq. ( |L9| ) is obviously one of 
those Is this coupling constant convenient for phenomenological descriptions using tree level 
diagrams in the infrared ? presumably no. If one aims at this phenomenology, as do the proponents 
of the "pinch technique" []l2|, 13], one could easily redefine a new eligible coupling constant by 



pulling a massive gluon propagator out for the gluon leg amputation of the ghost-ghost-gluon Green 
function used to build the coupling fll^. Thus a new (k 2 ) = ctT(k 2 )(k 2 +M 2 )/k 2 where M could be 



the gluon mass 2 . 



"We thank A.C. Aguilar, D. Binosi, J. Cornwall and J. Papavassiliou for this comment. 
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4. Can the bare ghost dressing function be finite non zero 



This question was raised by Kondo's remark [ 15, 16] of a relation between the k = values of 
the ghost dressing function F(k), Zwanziger's horizon function h(k), Kugo's function u(k) [ p7| , [03], 
and an additional function w{k). Applying to this relation Zwanziger's horizon gap equation and 
assuming that w(0) = he derives the surprising result that u(0) = —2/3 and F(0)=3, independently 
of the cut-off. The questions we will raise are: is this relation exact? does the prediction agree with 
lattice ? Is it compatible with renormalisability of QCD ? Our point of view is detailed in [jT^]. 

4.1 Kondo's relations 

In this subsection we only consider bare quantities. One solution to the problem of Gribov's 



ambiguity, which was proposed by Zwanziger ||20|], consists in using the Gribov-Zwanziger parti 



tion function, which aims at restricting the Gribov copies [21] within the Gribov Horizon: 



[DA] 8 (dA) det(M) e 



-s YM + yfd D xh(x) 



(4.1) 



for the D-dimensional Euclidean Yang-Mills theory, where SVm stands for the Yang-Mills action, 
M is the Faddeev-Popov operator, 



(4.2) 



M ab = -d^Df = (d^ + gf^A^ 
and h(x) is Zwanziger's horizon function, 

h(x) = J d D y gf hc Al{x){M- l r{x )y )gff e 4{y) ; 



(4.3) 



that restricts the integration over the gauge group to the first Gribov region, provided that the Gribov 
parameter, y, is a positive number. 

One defines then the function u(k 2 ) which, at vanishing momentum, gives the Kugo-Ojima 
parameter, and the function w(k 2 ) via the following identities. 



( (pfc b ) (gf de A d v c e ) ) k = - 8l v 8 ac u(k 2 ) ; 

(c a (gf def A e v cf) )l n = i8 ad k v (u(k 2 ) + w(k 2 )) . (4.4) 
From these definitions one obtains [FFl, Il6|] and rfl9l] without any hypothesis about u and w, 



b(0,A) 



F(0,A)-1 
D-l 



D 



D-l 



w(0,A) = -l-n(0,A) + 



1 



(MO))t=o 
D(A 2 -1) 

F(0,A) + (D-2) 



F(0,A) 



D-l 



+ 



1 



+ ■ 



D 



F(0,A) D-l 



(4.5) 

(h(0)) k=0 
D{N 2 -l) 



4.2 No finite F(0, A) is possible at large cut-off A 

If we use Zwanziger's gap equation: 



(A 2 -1) D 



(4.6) 
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4 5 6 

F(0,A) 



Figure 3: The solutions for u(0, A) and w(0, A) plotted as a function of F(0,A) under the assumption that 
the horizon gap equation is valid. 



K(A) 

u(0,A) 
• Kondo's solution 
co(0,A) 




5 6 

F(0,A) 



Figure 4: The same plot shown in fig. 3 but w(0, A) is required to be zero and the gap equation is relaxed by 
a multiplicative factor Jf(A), as explained in the text. fc(A) is plotted on the solid blue line. Again, current 
lattice estimates lie inside the green dotted square. 
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the functions «(0) and w(0) become, from eq. (4.5), a function of the bare F(0) plotted in Fig. [3]. 
The current lattice solutions for the bare ghost dressing functions at vanishing momentum lie in- 
side the green dotted square. The apparent approximate agreement of lattice results with Kondo 's 
solution is nevertheless misleading and due to the moderate cut-off value on the lattices. 

Indeed let us assume a fixed renormalised Fr(0, }X 2 ). The plot in fig. || can then be understood, 
as a function of the Zt,(pl 2 ,A) at fixed /I 2 as soon as inverse powers of A become negligible in front 



of logarithms, since, from eq. (|1.3[), eq. (1.4), 



Z 3 (At 2 ,A) 



*X0,A) 



+ 



,J_. 



(4.7) 



The large cut-off dependance of Z3 is known to be: 

,2 A\ / W/-A / A A \ 9/44 



Z 3 (M 2 ,Ao) Vlog(Ao/A Q c D )J 



(4.8) 



Z3(/i 2 ,A) — > 00 when A — > 00. 

Then the infinite cut-off limit is the limit at infinity on the horizontal axis of fig. |. The 
particular solution proposed in ref. [|l^, 16] (black circles), obtained by imposing w(0,A) = 0, 
corresponds to the intersection of u + w and u. It cannot hold when Z3 — > 00. Notice that the 
hypothesis of a finite bare F(0) with a vanishing F R (0) does not hold either since then F(0,A) = 
%F R {0) + <^(1/A 2 ) = G{\ /A 2 ) and consequently F(0, A) vanishes when A -> 00. 

Notice also from fig. || that a finite w(0, A) is not possible at the large cut-off limit. 

We should now take into account that gap equation (|4.6|) is a consequence of Gribov-Zwanziger 



modification of the Yang-Mills action eq. (4. 1 ). This is not what is done in LQCD, although lattice 



gauge fixing also restricts the Gribov copies within Gribov's horizon. Therefore we believe that 



condition eq. (4.6) has no reason to be fulfilled in LQCD and maybe not at all in QCD. Let us 
define Jc(A) such that 



(h(0)) k=0 = Jim^/ d D x (h{x))e ihx = k(A) (N 2 - l) D 



(4.9) 



If the gap equation eq. ( |4.6| ) is thus relaxed it becomes possible to keep w(0) finite, a result 
derived in [22, 23] in the Landau background gauge. We show the solution when w(0) =0 on 
the fig. |]. Nothing changes concerning the fact that the infinite cut-off limit is at infinity on the 
horizontal axis. Our conclusion still remains valid: is not possible to have a finite F(0,A) in the 
large A limit 

5. Conclusion 



In fig. |5] we perform an extrapolation of the lattice bare ghost dressing function using the 
small momentum expansion [ |19| ] shortly explained in eq. (2.4) of section 
larger lattice volumes are taken from ref. 



and the others from refs. [24, 



|. The data for the two 
. The fit with formula 



eq. (2.4) is rather good in its range of validity (small momentum). We also notice that the different 
lattice results seem to agree rather well although they correspond to different j8 values i.e. different 
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Figure 5: Bare ghost dressing function estimated from different lattice data sets. The solid line is for the 
best fit with the small-momentum expansion with R(P = 5.7(80 4 )) and the dashed one stands for the best fit 
withfl(/3 =5.7,oo). 

lattice spacings. One may feel happy and claim that we have a good scaling invariance. But this is 
a wrong statement from a misleading observation. 

Indeed, remember that the lattice spacing is the inverse of the cut-off in lattice regularisation. 
From eq. (<k8) we know that when the lattice spacing goes to zero, (A — > 00), F(0, A) °= /3 9//44 . On 



the whole range of lattice spacings considered in fig. |5| although the cut-off varies by more that a 
factor of 3, /3 9 / 44 varies only by 2.5%. This is why this variation is obscured by statistical errors 
in fig. |[ This fake "scaling invariance" hides the truth: F(k 2 ,A) rises very slowly to infinity when 
A — > °° i.e. /3 — ► °°. 

Bare values depend dramatically, although slowly, on the cut-off and have no real meaning 
unless the cut-off is specified. What makes really sense and has well defined limits at infinite 
cut-off (vanishing lattice spacing) are renormalised quantities [||, ||]. If we choose 1.5 GeV as the 
renormalisation scale, we get from lattice the gross estimate 

F (1.5 GeV) =Z 3 ~ 1.6 whence Fg (0, 1.5 GeV) ~ 2.2. (5.1) 

Altogether, combining all which has been discussed here, our conclusion concerning the ghost 
dressing function is 

• The renormalised ghost dressing function F#(0, /I 2 ) has a finite limit at vanishing momentum, 
F?(0, (1.5GeV) 2 ) ~ 2.2. It is a positive decreasing function at small momenta, probably also 
decreasing for all momenta. 



'indeed there is a trend of the largest cut-off, j3 = 6.4, to lie above the others, but this is hardly visible. 
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• The bare ghost dressing function F(k 2 ,A) goes very slowly to infinity at infinite A for all 
momenta. 
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